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(54) Title: BROADBAND ABSORPTION SPECTROSCOPY IN TURBID MEDIA BY COMBINED FREQUENCY-DOMAIN 
AND STEADY STATE METHODOLOGIES 

(57) Abstract: A technique for measuring broadband near-infrared absorption spectra of turbid media that uses a combination of 
frequency-domain and steady-state reflectance methods. Most of the wavelength coverage is provided by a white-light steady- state 
measurement, whereas the frequency -do main data are acquired at a few selected wavelengths. Coefficients of absorption and reduced 
scattering derived from the frequency -do main data are used to calibrate the intensity of the steady-state measurements and to deter- 
mine the reduced scattering coefficient at all wavelengths in the spectral window of interest. The absorption coefficient spectrum is 
determined by comparing the steady-state reflectance values with the predictions of diffusion theory, wavelength by wavelength. Ab- 
sorption spectra of a turbid phantom and of human breast tissue in vivo derived with the combined frequency -do main and steady-state 
technique, agree well with expected reference values. 
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BROADBAND ABSORPTION SPECTROSCOPY IN TURBID MEDIA BY 
COMBINED FREQUENCY-DOMAIN AND STEADY STATE METHODOLOGIES 

Background of the Invention 

Government Support 

The work described herein was supported by the National Institute of Health 
under grants GM50958 and RR01192 (Laser Microbeam and Medical Program), the 
U.S. Department of Energy (grant DE-FG03-91ER61227), and the U.S. Office of Naval 
Research (grant N00014-91-C-0134). 

1 . Field of the Invention 

This invention relates to methods and apparatus for analyzing turbid media by 
means of frequency-domain photon migration measurements and steady-state 
reflectance measurements. While not limited thereto, this invention is particularly useful 
for analyzing human body tissue in vivo. As one example, this invention is useful for 
analyzing female human breast tissue for detecting abnormal tissue conditions. 

2. Description of the Prior Art 

Reflectance spectroscopy is a technique for characterizing turbid media that has 
become widely used in medical diagnostics. In many cases quantification of 
chromophore concentrations is desired, and this requires the ability to separate the 
effects of absorption from those of scattering. Fundamentally, the coefficient of 
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absorption m a and the coefficient of reduced scattering p s ' can be determined by a 

series of reflectance measurements performed in one of three domains, namely, time 
(with a fast pulse of light), frequency (with a sinusoidally modulated source of light), and 
steady state (with a source of constant intensity but multiple detectors at different 
5 distances). 

Unsurprisingly, these three techniques have different merits and limitations. 
Spatially resolved steady-state techniques are relatively inexpensive and are more 
readily suited for the determination of p a and \i s ' over large, continuous ranges of 

wavelengths than are the other methods. However, the steady-state approach works 
10 best when measurements are performed with a combination of short (-1 transport mean 
free path) and long (many transport mean free paths) source-detector separations. 
Ideally, the optical properties of the sample should not vary over the ranges of volumes 
probed by the various measurements. The larger the spread of distances probed, the 
more likely that heterogeneities, such as those found in biological tissue, will distort the 
15 data from the predictions of the model. One approach to limiting this effect, given that 
the shortest separations provide great stability for the calculation of p s \ is to use 
relatively short (<10 mm) source-detector separations. Inasmuch as the mean 
probing depth scales with the source-detector separation, with this approach such 
measurements are sensitive to superficial components (to depths of less than 5 mm for 
20 typical biological tissues). 

Time-domain and frequency-domain techniques are well suited for deeper 
(>1 cm for biological tissue) investigations. Moreover, they can be performed with only 
one or a few source-detector separations, which makes them more robust for use in 
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studying heterogeneous samples. Because such techniques require sources that can 
be pulsed or modulated rapidly, covering a large wavelength range requires a tunable 
laser or an extensive collection of laser diodes, both of which can be expensive, difficult 
to maintain, and slow to cover the entire spectrum. This is an important drawback, 
because the quantification of chromophore concentrations can be significantly affected 
by use of a limited number of wavelengths. 

Brief Summary of the Invention 

The present invention provides a way to use steady-state and frequency-domain 
reflectance measurements in tandem to obtain broad wavelength coverage with 
increased penetration depth. In accordance with one feature of the invention, there is 
provided a method of analyzing turbid media. This method comprises performing 
frequency-domain photon migration measurements on a turbid medium sample. It also 
includes performing steady-state reflectance measurements on the same turbid medium 
sample. It further includes combining the results of the frequency-domain photon 
migration measurements with the results of the steady-state reflectance measurements 
for obtaining a unique spectrum for the turbid medium sample. 

In accordance with another feature of the invention, there is provided a method of 
analyzing human body tissue in vivo. This method includes emitting light photons from 
a plurality of different laser diode light sources into human body tissue from a 
predetermined delivery point on the surface of a human body. This method further 
includes emitting light photons from a white light source into the same human body 
tissue from the same predetermined delivery point on the surface of the human body. 
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This method also includes collecting the light photons received at a spaced collection 
point on the surface of the human body after such fight photons have traveled through 
the human body tissue intermediate the delivery and collection points. This method 
further includes combining predetermined characteristics of the light photons from the 
different light sources to provide an indication of the composition of the traversed 
human body tissue. . 

While the apparatus and method has or will be described for the sake of 
grammatical fluidity with functional explanations, it is to be expressly understood that 
the claims, unless expressly formulated under 35 USC 112, are not to be construed as 
necessarily limited in any way by the construction of "means" or "steps" limitations, but 
are to be accorded the full scope of the meaning and equivalents of the definition 
provided by the claims under the judicial doctrine of equivalents, and in the case where 
the claims are expressly formulated under 35 USC 112 are to be accorded full statutory 
equivalents under 35 USC 112. 

For a better understanding of the present invention, together with other and 
further advantages and features thereof, reference is made to the following description 
taken in connection with the accompanying drawings, the scope of the invention being 
pointed out in the appended claims. 

Brief Description of the Drawings 

FIG. 1 is a schematic diagram of a representative embodiment of apparatus for 
practicing the present invention; 

FIG. 2 is a chart used for explaining a representative embodiment of the method 
of the present invention; 
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FIG. 3 is a graph of a steady-state reflectance spectrum obtained for a dye- 
Intralipid phantom; 

FIGS. 4-6 are graphs of steady-state reflectance spectra acquired from three 
locations in breast tissue of female volunteers; 

FIG. 7 is a graph of a broadband reduced scattering coefficient spectrum 
determined for the same dye-Intralipid phantom used for FIG. 3; 

FIGS. 8-10 are graphs of broadband reduced scattering coefficient spectra 
determined for the three female breast locations of FIGS. 4-6, respectively; 

FIG. 11 is a graph of a broadband absorption coefficient spectrum determined for 
the same dye-Intralipid phantom used for FIG. 3; 

FIGS. 12-14 are graphs of broadband absorption coefficient spectra determined 
for the three female breast locations of FIGS. 4-6, respectively; and 

FIGS. 15-18 are graphs of absorption coefficient spectra of major absorbers in 
breast tissue. 

The invention and its various embodiments can now be better understood by 
turning to the following detailed description of the preferred embodiments which are 
presented as illustrated examples of the invention defined in the claims. It is expressly 
understood that the invention as defined by the claims may be broader than the 
illustrated embodiments described below. 

Detailed Description of the Preferred Embodiments 

Referring to FIG. 1, there is shown a schematic diagram of a representative 
embodiment of apparatus 10 constructed in accordance with the present invention for 
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analyzing turbid media. Such apparatus 10 includes apparatus 11 for performing 
frequency-domain photon migration (FDPM) measurements on a turbid medium sample 
12. The frequency-domain (FD) apparatus 1 1 includes a plurality of laser diodes 13 for 
emitting light energy at a plurality of different wavelengths. In the present embodiment, 

5 the plurality of laser diodes 1 3 includes seven laser diodes which emit light energy at 
wavelengths of 672, 800, 806, 852, 896, 913 and 978 nanometers (nm). At least one 
delivery optical fiber 14 is provided for delivering modulated light energy from the laser 
diodes 13 to a delivery point on the surface of the turbid medium sample 12. The output 
powers of each of laser diodes 13 is less than 20 milliwatts (mW) at the sample 12. 

10 A mechanism (not shown) is provided for modulating the light energy emitted by 

the laser diodes 13. This modulating mechanism may take the form described in U.S. 
patent 5,424,843, granted on June 13, 1995 to Bruce J. Tromberg et al. This patent 
describes a form of apparatus for making frequency-domain measurements and the 
descriptions set forth in this patent are hereby incorporated herein by this reference 

15 thereto. 

The frequency-domain apparatus 1 1 further includes a collection optical fiber 15 
for collecting emitted light energy from the laser diodes 13 arriving at a collection point 
on the surface of the turbid medium sample 12 after traveling through the sample 12. 
The collection point is fixed or located at a predetermined distance "r" from the delivery 
20 point. This distance r is, for example, in the range of 15 to 25 millimeters (mm). 

The frequency-domain apparatus 1 1 also includes an analyzer mechanism 
coupled to the collection optical fiber 15 for measuring at least one characteristic of the 
modulated light energy received from the laser diodes 13 after passage through the 
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turbid medium sample 12. This mechanism includes an avalanche photodiode (APD) 
unit 16 that detects the light photons collected by the collection optical fiber 15 and 
amplifies the alternating-current component of the detected signals. A network analyzer 
(e.g., a Hewlett-Packard model 8753C) is coupled to the avalanche photodiode unit 16 

5 and measures phase and modulation amplitudes of the detected photon signals. 

The analyzing apparatus 10 further includes apparatus 20 for performing steady- 
state (SS) reflectance measurements on the same turbid medium sample 12. This 
steady-state apparatus 20 includes a white light source 21 for emitting white light 
covering the visible spectrum from at least 500 to at least 1200 nm. A suitable white 

10 light source is a 150-watt halogen electric lamp. Light energy from the white light 

source 21 is delivered to the turbid medium sample 12 by way of the delivery optical 
fiber 14. 

The steady-state (SS) apparatus 20 also includes a spectrograph mechanism 22 
coupled to the collection optical fiber 15 for receiving white light energy arriving at the 

15 collection point on the turbid medium sample 12 and producing therefrom a reflectance 
spectrum for the sample 12. Spectrograph mechanism 22 may be, for example, an 
Ocean Optics model S2000 with a linear CCD detector from 525 to 1 1 55 nm. This 
Ocean Optics spectrograph records a total of 2048 points (0.35 nm/pixel) with a spectral 
resolution of 5 nm. Relative reflectance is calculated to be the sample value divided by 

20 the source value. 

When making measurements, the various light sources are operated in a 
sequential manner. Light energy is emitted by the different laser diodes, one laser 
diode at a time. Then white light alone is emitted by the white light source for a short 
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time interval. Total measurement time per sample is on the order of 40 seconds (30 
seconds for the laser diode measurements and 10 seconds for the white light 
measurements). 

The measured reflectance signal R can be predicted by use of the diffusion 
approximation to the radiative transfer equation- In this theoretical framework, the 
reflectance is a function of the optical properties of the medium, defined by the 
absorption coefficient p a , the reduced scattering coefficient p s \ and the index of 

refraction "n". The diffusion approximation is valid for large source-detector separation r 

[r > 10(p~ + Ms')" 1 ] and h 'g h (reduced) albedo [ps7 (p + p s ') > 0.95]. The steady-state 

and frequency-domain cases can be described by a single formalism in which the 
solution for the reflectance is built from the Green's function for the diffusion equation, 

i.e., the fluence (W/cm 2 ) that is due to an isotropic point source in an infinite, 
homogeneous medium. This function takes the form exp(-kr)/(Dr), where D = [3(p a + 

Ms')]" 1 , the complex wave number is k = /c rea | + /^i ma ge ■ ancl 



CO 



+1 



1 + 
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where 00 is the modulation frequency in radians per second, and c is the speed of light. 
The steady-state solution is simply the limit when w = 0. In this case, the solution for R 

is a real number {/c rea [ = [3p a (M a + Ms')]^ ancl ^image = °}- When w > 0, R becomes a 
complex number A[exp(~/0)], where A is the modulation amplitude and O is the 
phase shift relative to the source. These quantities are obtained from the frequency- 
domain measurement. 

In treating reflectance problems, the sample-air interface is modeled using an 
extrapolated boundary condition in which the fluence is set to zero at a distance z^ = 

2D(1 + R e ff )/(1 - R e ff ) above the sample R e ff depends on the refractive mismatch and 
equals 0.493 for tissue of n = 1 .4 and air of n = 1.0). The method of images is 
employed, with an isotropic point source at a depth z Q - ^ + p s ')-1 contributing a 

signal S r and a negative image point at a height z Q + z^ above the extrapolated 

boundary contributing Sj . The result is the fluence O at any point in the sample. The 

detected signal along the boundary is then written as a combination of terms 
proportional to the fluence and to its flux normal to the surface: 

where values for the constants c 1 and c 2 are determined by the refractive index 

mismatch between the two media (for tissue of n = 1 .40 and air of n = 1 , for these 
constants the values 0.1 18 and 0.306, respectively, are assumed), -z is a unit vector 
pointing normally upward out of the sample, and 
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exp(-kr s ) exp(-&r) 
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where P is the incident power and r s is the distance from the source (image) to the 

detector- 
Each frequency-domain measurement made with each laser diode contains 
instrumental artifacts. The measured reflectance is actually 

R = C Q A exp[- / (cp + cp 0 M ( 3 ) 

where C 0 and cp 0 are sample-independent instrumental constants. A and cp are, 

respectively, the modulation amplitude and phase of the frequency-domain reflectance, 
as defined above. 

For calibration, frequency-domain data is gathered from a prepared sample 
whose p a and p s ' values are known a priori from a set of two-distance frequency- 
domain measurements. Frequency-dependent values of cp 0 and C Q are calculated 

from the discrepancies between measured [Equation (3)] and predicted [Equation (1)] 
phase and modulation amplitude, thus calibrating future measurements. 

With the instrumental constants thus determined, the 502 data points per sample 
(251 for both phase and amplitude as a function of the frequency-domain modulation 
frequency w) depend on two unknowns, p a and (j s '. Selection is made of the best p a 
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and |j s ' values to fit the predictions of Equation (1) to the data, using the iterative, non- 
linear, least-squares method of Levenberg and Marquardt , simultaneously fitting the 
phase and the amplitude until overall convergence is achieved- The reflectance in 
Equation (1) is nearly proportional to either the flux or the fluence term alone (i.e., the 
ratio of the two is independent of p a and p s ') for the large source-detector separations 

employed in the illustrated embodiment [>10(|J a + Ms')" 1 ]- Consequently, only the 

fluence term O from Equation (1) is used herein, the additional proportionality factor 
being absorbed into the instrumental constant C 0 . 

Additional information on the foregoing theory is set forth in a paper by T.H. 
Pham, O. Coquoz, J.B. Fishkin, E. Anderson and B.J. Tromberg entitled "Broad 
bandwidth frequency domain instrument for quantitative tissue optical spectroscopy", 
Rev. Sci. Instrum., vol. 71, pages 2500-2513, (2000). The contents of this paper are 
hereby incorporated herein by this reference thereto. Consideration will now be 

given in regard to how frequency-domain and steady-state measurements can be 
combined to yield quantitative broadband absorption coefficient (p a ) spectra. The goal 

is to compute p a at each wavelength, given the measured steady-state reflectance. 

However, single-distance steady-state reflectance cannot itself provide p a . R depends 

on Ms' as well as on pa, so one measurement cannot provide a unique determination of 
either parameter. In addition, the instrumental constant C Q is not known. The single 
steady-state measurement at each wavelength therefore needs to be supplemented 
with both a means of establishing the absolute reflectance intensity and an additional 
piece of information about p a or p s ' (or a combination of the two). 
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The frequency-domain system can provide both pieces of necessary information 
at each steady-state wavelength over the continuous steady-state spectrum. At first 
glance, this is surprising because the frequency-domain system operates at only seven 
discrete wavelengths, whereas the steady-state system covers a 450 nm range with 
continuous 5-nm resolution. It would seem that the needed values of \j s ' and C Q could 

be determined only at those these seven wavelengths, and strictly speaking such is 
indeed the case. Significantly, however, the wavelength dependence of both jj s ' and C 

0 is smooth and predictable in shape. It is this fact that enables use of the discrete 

frequency-domain information to supplement the broadband steady state 
measurements at all desired wavelengths. 

The frequency-domain measurements permit calculation of the instrumental 
factor C Q at all wavelengths because no wavelength dependence at ail is expected. 

The sample and source spectra are measured with the same delivery and collection 
system and only a few seconds apart, so there should be negligible wavelength- 
dependent artifacts in the ratio of the two. Therefore the task reduces to that of 
calculating C Q at a single wavelength. This is readily done at any of the frequency- 
domain wavelengths, as Ma and Ms' are both known and absolute reflectance is a 
function only of these variables [cf. Equations (1) and (2) in the co = 0 limit]. Using all 
the frequency-domain data to increase robustness, we can calculate the value of C Q 
that scales the measured steady-state reflectances to match the predicted reflectances 
as closely as possible (in the least-squares sense). Once this is done, the scale factor 
C G is known for the entire spectrum. 
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The frequency-domain values of jj s ' also allow us to obtain information about p s ' 
across the entire wavelength range. The particle size distribution of scatterers (0.1-10 
mm) in many biological media and phantoms tend to have smooth wavelength 
dependence over the wavelength range 650-1000 nm which is well described by a 
power function of the form: 

Ms '(A)=AA~B (4) 
[the parameter A in Equation (4) is a constant and has no relation to the modulation 
amplitude]. As a result, measuring several p s ' values allows us to fit those values to a 
simple function of wavelength and to obtain good estimates of p s ' at all other 
wavelengths needed. 

Referring to FIG. 2, there is shown a flowchart for the overall process. As 
indicated by block 24, the frequency-domain photon migration (FDPM) measurements 
on the sample provide determinations of the absorption coefficient (p a ) values and the 

reduced scattering coefficient (n s ' ) values at the discrete wavelengths A ^- A 5 of the 

light energy emitted by the different ones of laser diodes 13. As indicated by block 25, 
the steady-state (SS) reflectance measurements on the sample provide a determination 
of a relative reflectance spectrum aR(A) for the sample for the light energy emitted by 
the white light source 21 , where a denotes a proportionality constant. For the present 
example, the white light energy and, hence, the reflectance spectrum is assumed to 
cover a wavelength range of 650-1000 nm. 

First, the reflectance R predicted by steady-state diffusion theory is calculated 
at each discrete one of frequency-domain wavelengths A 1 - A 6 , using the p a and p s ' 
values supplied by the FDPM measurements 24. This is indicated at block 26. These 
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six reflectance values are absolute reflectance values. By scaling the relative 
reflectance curve provided by steady-state measurements 25 to fit these six absolute 
values in a least squares sense, the proportionality factor a is removed and the 
reflectance spectrum curve is converted to absolute units. This is indicated by block 27, 
where R(A) denotes the absolute reflectance spectrum. 

Next, the scattering coefficient spectrum p s ' (A) is deduced from the six p s ' 
frequency-domain measurements. This is indicated by block 28. This is possible 
because the particle size distribution of scatterers in the samples studied tends to have 
smooth wavelength dependence over the 600-1000 nm range which is well described 
by a power function of the form set forth in Equation (4). Consequently, the six 
measured scattering coefficient values are fitted to the matching power law curve to 
provide a complete scattering coefficient spectrum curve for the 600-1000 nm range. 

As indicated by block 29, the results of the frequency-domain photon migration 
measurements are combined with the results of the steady-state reflectance 
measurements to obtain an absorption coefficient spectrum p a (A) for the entire 600- 

1000 nm range. More particularly, the absorption coefficient spectrum value is 
numerically calculated at each wavelength from the corresponding calibrated or 
absolute reflectance value obtained from the calibrated reflectance spectrum curve and 
the corresponding scattering coefficient value obtained from the scattering coefficient 
spectrum curve. The p a (A) values were calculated with an in-house Matlab (The Math 
Works, Inc.) program, making use of the optimization toolbox. Solving for the entire p a 

spectrum required approximately ten seconds on a personal computer. 
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The above-described method of combining the steady-state and frequency- 
domain methods was tested by using it to measure the absorption spectra of turbid 
samples- Analysis was made of a liquid tissue phantom whose absorption spectrum is 
known by direct spectrophotometry of the absorbing component before mixing. In vivo 
measurements were also made on the breasts of two human female volunteers, 
demonstrating that data can be gathered in vivo and analyzed to provide pertinent 
physiological parameters. 

The phantom contained 17 mg of green dye (naphthol) and 50 mL of an aqueous 
scattering suspension (lntralipid-20%, Pharmacia, Inc.) dissolved in 950 mL of water. 
Before the addition of Intralipid, the absorbance spectrum of the dye-water solution was 
measured in a 1-cm plastic cuvette by a spectrophotometer (Beckman Instruments, 
DU630) with water as a reference. Optical measurements were performed as indicated 
in FIG. 1 , with the source and detector optical fibers placed at the surface of the sample. 

In vivo measurements were performed on the left breasts of two supine female 
volunteers, aged 37 and 21 years. Data were gathered from two regions on one 
volunteer (a region on the areolar border and a region close to the middle of the chest) 
and from one region (inner breast only) on the other volunteer, yielding a total of three 
samples. The source light was delivered by optical fiber, but the frequency-domain 
detector was placed directly against the tissue, without a collection fiber. The source 
fiber and detector were bundled into a single hand-held probe that was placed gently 
against the breast. The steady-state reflectance was measured subsequently, at the 
exact same location, in the two-fiber mode shown in FIG. 1. 
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FIG. 3 is a graph showing the steady-state reflectance measured from the dye- 
Intralipid phantom (solid curve) and the predicted absolute reflectance calculated from 
diffusion theory based on the discrete frequency-domain measurements of M a anc l Ms' 

(solid dots). Because of strong absorption at wavelengths above 950 nm, the 
frequency-domain measurement at 978 nm was noisy, which in turn made the 
calculation of steady-state data at wavelengths above 950 nm unstable. Data is 
therefore presented for 650-950 nm for the phantom only. As anticipated, all the 
measurements differ from the predictions by essentially the same scale factor. The 
entire steady-state spectrum is thus converted into absolute reflectance units. 

Similar results were obtained for the three female breast samples, as shown in 
FIGS. 4-6. FIG. 4 is from the inner breast region of the 37-year-old female and FIGS. 5 
and 6 are from the areolar border and the inner breast, respectively, of the 21 -year-old 
female. 

FIG. 7 shows the steady-state scattering coefficient spectrum (solid curve) 
determined for the dye-lntralipid phantom. Open circles represent the discrete values 
measured by the frequency-domain technique and the solid curve shows the best 
power-law fit to Equation (4). Fitting parameters are A = 2200 and B = - 0.82. Error 
bars on the discrete values are 3% bars. As expected, the wavelength dependence is 
fairly smooth and is easily described by the fit. Corresponding fits for the three female 
breast samples are shown in FIGS. 8-10. Fitting parameters (A,B) are, respectively, 
(240, - 0.86), (2700, - 12) and (250, - 0.85). Note that the values change by only a 
factor of two over the entire range, so the wavelength dependence is weak in addition to 
being predictable. 
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FIG. 1 1 is a graph showing the absorption coefficient (m a ) spectrum calculated 

for the dye-lntralipid phantom- The thicker curve shows the absorption coefficient 
values that were found to make Equation (1) reproduce the measured reflectance data 
as closely as possible, using the zero algorithm as described above. The thinner curve 
represents a linear least-squares fit to the thicker curve based on reference spectra of 
pure water and pure dye (measured by a spectrophotometer). As FIG. 1 1 shows, the 
two-component fit accurately models the measured spectrum across the entire spectral 
range. Also shown are the values of p a derived by frequency-domain analysis alone 

(open circles), demonstrating an agreement between the frequency-domain and the 
combined FD-SS methods at these wavelengths. 

Corresponding plots that show absorption coefficient spectra for the three female 
breast samples appear in FIGS. 12-14. As before, the thicker curve is the experimental 
data, the thinner curve is a full spectrum fit and the open circles are frequency-domain 
values. In these cases, the fit is built from published spectra (FIGS. 15-18) of 
oxyhemoglobin, deoxyhemoglobin, water and fat, which are commonly regarded as the 
four major absorbers in breast tissue in the 650-1000 nm range. Coefficients from the 
fit thus provide estimates of these component concentrations. 

For comparison purposes, the dotted curves in FIGS. 12-14 show the spectral fit 
obtained when only the frequency-domain wavelengths are weighted, as in a customary 
FD-only experiment. The two spectral reconstructions are clearly different, and these 
different reconstructions lead to different concentration estimates, 

FIG. 12 is for the inner breast measurement on the 37-year-old female. 
Physiological parameters from this spectral fit are: total hemoglobin concentration, 
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22 mM; oxygen saturation, 73%; water, 15 g/cm^; fat, 0-75 gm/cm 0 *. 

FIG, 13 is for the areolar border breast region of the 21 -year-old female. 
Physiological parameters from FIG. 13 are: total hemoglobin concentration, 30 mM ; 

oxygen saturation, 70%; water, 0.51 g/cm 3 ; fat, 0.42 g/cm 3 . 

FIG. 14 is for the inner breast of the 21 -year-old female. Physiological 
parameters for FIG. 14 are: total hemoglobin concentration, 19 mM; oxygen saturation, 

72%; water, 28 g/cm 3 ; fat, 0.56 g/cm 3 . 

FIGS. 15-18 show the absorption coefficient spectra of major absorbers in breast 
tissue. FIG. 15 shows the spectrum for oxyhemoglobin, 1 mM. FIG. 16 shows the 

spectrum for deoxyhemoglobin, 1 mM . FIG. 17 shows the spectrum for water, 1 g/cm 3 . 

FIG. 18 shows the spectrum for fat (soybean oil), 0.9 g/cm 3 . 

Compared with discrete wavelength frequency-domain measurements alone, 
access to a continuous absorption coefficient spectrum offers two important 
advantages, namely, better chromophore identification and improved concentration 
quantification. The female breast spectra (FIGS. 12-14) illustrate the usefulness of the 
combined frequency-domain and steady-state measurements for clinical investigations. 
Among other things, the resulting indications of tissue composition can be used to 
detect abnormal tissue conditions. 

In summary, the invention can thus be defined as a method of analyzing turbid 
media comprising the steps of: performing frequency-domain photon migration 
measurements on a turbid medium sample; performing steady-state reflectance 
measurements on the same turbid medium sample; and combining the results of the 
frequency-domain photon migration measurements with the results of the steady-state 
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reflectance measurements for obtaining a unique spectrum for the turbid medium 
sample. 

The unique spectrum is a calibrated reflectance spectrum for the turbid medium 
sample. 

The steady-state reflectance measurements provide a relative reflectance 
spectrum and the calibrated reflectance spectrum is obtained by scaling the relative 
reflectance spectrum to fit absolute reflectance values determined at a plurality of 
discrete wavelengths by the frequency-domain photon migration measurements. 

A photon scattering coefficient spectrum is determined from the frequency- 
domain photon migration measurements, the reflectance measurements provide a 
reflectance spectrum and the absorption coefficient spectrum is calculated from the 
scattering coefficient spectrum and the reflectance spectrum. 

The frequency-domain photon migration measurements are made at a plurality of 
discrete wavelengths. The frequency-domain photon migration measurements are 
made by sequentially exposing the sample to light emitted by a plurality of laser diodes 
which emit light at different wavelengths. The steady-state reflectance measurements 
are made by exposing the sample to white light emitted by a white light source. 

The invention is also characterized as a method of analyzing turbid media 
comprising the steps of: performing frequency-domain photon migration measurements 
on a turbid medium sample at a plurality of discrete wavelengths for determining photon 
absorption coefficient values and photon scattering coefficient values of the turbid 
medium at such discrete wavelengths; performing steady-state reflectance 
measurements on the same turbid medium sample for determining a relative reflectance 
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spectrum for the sample; calibrating the relative reflectance spectrum obtained from the 
steady-state measurements in accordance with the discrete absorption and scattering 
coefficient values obtained from the frequency-domain measurements; determining a 
scattering coefficient spectrum from the discrete scattering coefficient values obtained 
from the frequency-domain measurements; and combining calibrated reflectance 
spectrum values and scattering coefficient spectrum values to obtain an absorption 
coefficient spectrum for the turbid medium sample- 
Again the relative reflectance spectrum is calibrated by calculating from the 
frequency-domain photon migration measurements absolute reflectance values at the 
discrete wavelengths and scaling the relative reflectance spectrum to fit these absolute 
reflectance values. 

The scattering coefficient spectrum is determined by fitting the scattering 
coefficient values determined at the discrete wavelengths to a matching power law 
curve for the desired wavelength range. 

The absorption coefficient spectrum is determined by numerical calculation at 
each wavelength of the absorption coefficient value based on the calibrated reflectance 
value and the scattering coefficient value obtained from the calibrated reflectance 
spectrum and scattering coefficient spectrum curves. 

The invention can be described as a method of analyzing human body tissue in 
vivo comprising: emitting light photons from a plurality of different laser diode light 
sources into human body tissue from a predetermined delivery point on the surface of a 
human body; emitting light photons from a white light source into the same human body 
tissue from the same predetermined delivery point on the surface of the human body; 
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collecting the light photons received at a spaced collection point on the surface of the 
human body after such light photons have traveled through the human body tissue 
intermediate the delivery and collection points; and combining predetermined 
characteristics of the collected light photons from the different light sources to provide 
an indication of the composition of the traversed human body tissue. 

The indication of the composition of the traversed human body tissue takes the 
form of a photon absorption coefficient spectrum for the traversed human body tissue. 
The light photons from the different light sources are applied in a sequential manner. 
The human body tissue is female human breast tissue, the delivery and collection points 
are located on the surface of a female human breast and the resulting indication of 
tissue composition can be used to detect an abnormal tissue condition. 

The invention is still further defined as an apparatus for analyzing turbid media 
comprising apparatus for performing frequency-domain photon migration measurements 
on a turbid medium sample; apparatus for performing steady-state reflectance 
measurements on the same turbid media sample; and apparatus for combining results 
of the frequency-domain photon migration measurements with results of the steady- 
state reflectance measurements for producing an absorption coefficient spectrum for the 
turbid medium sample. 

The apparatus for performing frequency-domain photon migration measurements 
includes: a plurality of laser diodes for emitting light energy at a plurality of different 
wavelengths; a delivery optical fiber for delivering light energy from the laser diodes to 
the turbid medium sample; and a mechanism for receiving light energy arriving at a 
collection point on the turbid medium sample and measuring a characteristic of the 
collected light energy. 
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The apparatus for performing steady-state measurements includes: a white light 
source for emitting white light; a delivery optical fiber for delivering light energy from the 
white light source to the turbid medium sample; and a spectrograph mechanism for 
receiving white light energy arriving at a spaced collection point on the turbid medium 
sample and producing therefrom a reflectance spectrum for the turbid medium sample. 

The invention is still further summarized as an apparatus for analyzing turbid 
media comprising: a plurality of laser diodes for emitting light energy at a plurality of 
different wavelengths; a mechanism for modulating the light energy emitted by the laser 
diodes; a white light source for emitting white light; at least one delivery optical fiber for 
delivering modulated light energy from the laser diodes and white light from the white 
light source to a delivery point on a surface of a turbid medium sample; a collection 
optical fiber for collecting emitted light energy arriving at a collection point on the 
surface of the turbid medium sample after traveling through the turbid medium sample, 
such collection point being fixed at a predetermined distance from the delivery point; an 
analyzer mechanism coupled to the collection optical fiber for measuring at least one 
characteristic of the modulated light energy received from the laser diodes after 
passage through the turbid medium sample; a spectrograph mechanism coupled to the 
collection optical fiber for producing a reflectance spectrum for the white light received 
from the white light source after passage through the turbid medium sample; and a data 
processing mechanism for combining results of the analyzer mechanism measurements 
with results of the spectrograph mechanism measurements for producing a photon 
absorption coefficient spectrum for the turbid medium sample. 

Many alterations and modifications may be made by those having ordinary skill in 
the art without departing from the spirit and scope of the invention. Therefore, it must 
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be understood that the illustrated embodiment has been set forth only for the purposes 
of example and that it should not be taken as limiting the invention as defined by the 
following claims. For example, notwithstanding the fact that the elements of a claim are 
set forth below in a certain combination, it must be expressly understood that the 
invention includes other combinations of fewer, more or different elements, which are 
disclosed in above even when not initially claimed in such combinations. 

The words used in this specification to describe the invention and its various 
embodiments are to be understood not only in the sense of their commonly defined 
meanings, but to include by special definition in this specification structure, material or 
acts beyond the scope of the commonly defined meanings. Thus if an element can be 
understood in the context of this specification as including more than one meaning, then 
its use in a claim must be understood as being generic to all possible meanings 
supported by the specification and by the word itself. 

The definitions of the words or elements of the following claims are, therefore, 
defined in this specification to include not only the combination of elements which are 
literally set forth, but all equivalent structure, material or acts for performing substantially 
the same function in substantially the same way to obtain substantially the same result. 
In this sense it is therefore contemplated that an equivalent substitution of two or more 
elements may be made for any one of the elements in the claims below or that a single 
element may be substituted for two or more elements in a claim. Although elements 
may be described above as acting in certain combinations and even initially claimed as 
such, it is to be expressly understood that one or more elements from a claimed 
combination can in some cases be excised from the combination and that the claimed 
combination may be directed to a subcombination or variation of a subcombination. 

Insubstantial changes from the claimed subject matter as viewed by a person 
with ordinary skill in the art, now known or later devised, are expressly contemplated as 
being equivalent^ within the scope of the claims. Therefore, obvious substitutions now 
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or later known to one with ordinary skill in the art are defined to be within the scope of 
the defined elements. 

The claims are thus to be understood to include what is specifically illustrated 
and described above, what is conceptionally equivalent, what can be obviously 
5 substituted and also what essentially incorporates the essential idea of the invention. 
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We claim: 



1 1 . A method of analyzing turbid media comprising: 

2 performing frequency-domain photon migration measurements on a turbid 

3 medium sample; 

4 performing steady-state reflectance measurements on the same turbid medium 

5 sample; 

6 and combining the results of the frequency-domain photon migration 

7 measurements with the results of the steady-state reflectance measurements for 

8 obtaining a unique spectrum for the turbid medium sample. 

1 2. A method in accordance with Claim 1 wherein the unique spectrum is a 

2 calibrated reflectance spectrum for the turbid medium sample. 

1 3. A method in accordance with Claim 2 wherein the steady-state reflectance 

2 measurements provide a relative reflectance spectrum and the calibrated reflectance 

3 spectrum is obtained by scaling the relative reflectance spectrum to fit absolute 

4 reflectance values determined at a plurality of discrete wavelengths by the frequency- 

5 domain photon migration measurements. 

l 4. A method in accordance with Claim 1 wherein the unique spectrum is a 
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2 photon absorption coefficient spectrum for the turbid medium sample. 

1 5. A method in accordance with Claim 4 wherein a photon scattering coefficient 

2 spectrum is determined from the frequency-domain photon migration measurements, 

3 the reflectance measurements provide a reflectance spectrum and the absorption 

4 coefficient spectrum is calculated from the scattering coefficient spectrum and the 

5 reflectance spectrum. 

1 6. A method in accordance with Claim 1 wherein the frequency-domain photon 

2 migration measurements are made at a plurality of discrete wavelengths. 

1 7. A method in accordance with Claim 1 wherein the frequency-domain photon 

2 migration measurements are made by sequentially exposing the sample to light emitted 

3 by a plurality of laser diodes which emit light at different wavelengths. 

1 8. A method in accordance with Claim 1 wherein the steady-state reflectance 

2 measurements are made by exposing the sample to white light emitted by a white light 

3 source. 

1 9. A method of analyzing turbid media comprising: 

2 performing frequency-domain photon migration measurements on a turbid 

3 medium sample at a plurality of discrete wavelengths for determining photon absorption 
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4 coefficient values and photon scattering coefficient values of the turbid medium at such 

5 discrete wavelengths; 

6 performing steady-state reflectance measurements on the same turbid medium 

7 sample for determining a relative reflectance spectrum for the sample; 

8 calibrating the relative reflectance spectrum obtained from the steady-state 

9 measurements in accordance with the discrete absorption and scattering coefficient 
10 values obtained from the frequency-domain measurements; 

n determining a scattering coefficient spectrum from the discrete scattering 

12 coefficient values obtained from the frequency-domain measurements; 

13 and combining calibrated reflectance spectrum values and scattering coefficient 

14 spectrum values to obtain an absorption coefficient spectrum for the turbid medium 

15 sample. 

1 10. A method in accordance with Claim 9 wherein the relative reflectance 

2 spectrum is calibrated by calculating from the frequency-domain photon migration 

3 measurements absolute reflectance values at the discrete wavelengths and scaling the 

4 relative reflectance spectrum to fit these absolute reflectance values. 

1 1 1 . A method in accordance with Claim 9 wherein the scattering coefficient 

2 spectrum is determined by fitting the scattering coefficient values determined at the 

3 discrete wavelengths to a matching power law curve for the desired wavelength range. 

l 12. A method in accordance with Claim 9 wherein the absorption 
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2 coefficient spectrum is determined by numerical calculation at each wavelength of the 

3 absorption coefficient value based on the calibrated reflectance value and the scattering 

4 coefficient value obtained from the calibrated reflectance spectrum and scattering 

5 coefficient spectrum curves. 



1 1 3- Apparatus for analyzing turbid media comprising: 

2 apparatus for performing frequency-domain photon migration measurements on 

3 a turbid medium sample; 

4 apparatus for performing steady-state reflectance measurements on the same 

5 turbid media sample; 

6 and apparatus for combining results of the frequency-domain photon migration 

7 measurements with results of the steady-state reflectance measurements for producing 

8 an absorption coefficient spectrum for the turbid medium sample. 



1 14. Apparatus in accordance with Claim 13 wherein the apparatus for 

2 performing frequency-domain photon migration measurements includes: 

i 3 a plurality of laser diodes for emitting light energy at a plurality of different 

4 wavelengths; 

5 a delivery optical fiber for delivering light energy from the laser diodes to the 

6 turbid medium sample; 
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7 and a mechanism for receiving light energy arriving at a collection point on the 

8 turbid medium sample and measuring a characteristic of the collected light energy. 

1 1 5. Apparatus in accordance with Claim 1 3 wherein the apparatus for 

2 performing steady-state measurements includes: 

3 a white light source for emitting white light; 

4 a delivery optical fiber for delivering light energy from the white light source to the 

5 turbid medium sample; 

6 and a spectrograph mechanism for receiving white light energy arriving at a 

7 spaced collection point on the turbid medium sample and producing therefrom a 

8 reflectance spectrum for the turbid medium sample. 

1 16. Apparatus for analyzing turbid media comprising: 

2 a plurality of laser diodes for emitting light energy at a plurality of different 

3 wavelengths; 

4 a mechanism for modulating the light energy emitted by the laser diodes; 

5 a white light source for emitting white light; 

6 at least one delivery optical fiber for delivering modulated light energy from the 

7 laser diodes and white light from the white light source to a delivery point on a surface 

8 of a turbid medium sample; 

9 a collection optical fiber for collecting emitted light energy arriving at a collection 
10 point on the surface of the turbid medium sample after traveling through the turbid 

n medium sample, such collection point being fixed at a predetermined distance from the 
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12 delivery point; 

13 an analyzer mechanism coupled to the collection optical fiber for measuring at 

14 least one characteristic of the modulated light energy received from the laser diodes 

15 after passage through the turbid medium sample; 

16 a spectrograph mechanism coupled to the collection optical fiber for producing a 

17 reflectance spectrum for the white light received from the white light source after 
is passage through the turbid medium sample; 

19 and a data processing mechanism for combining results of the analyzer 

20 mechanism measurements with results of the spectrograph mechanism measurements 

21 for producing a photon absorption coefficient spectrum for the turbid medium sample. 

1 17. A method of analyzing human body tissue in vivo comprising: 

2 emitting light photons from a plurality of different laser diode light sources into 

3 human body tissue from a predetermined delivery point on the surface of a human 

4 body; 

5 emitting light photons from a white light source into the same human body tissue 

6 from the same predetermined delivery point on the surface of the human body; 

7 collecting the light photons received at a spaced collection point on the surface of 

8 the human body after such light photons have traveled through the human body tissue 

9 intermediate the delivery and collection points; 

10 and combining predetermined characteristics of the collected light photons from 

n the different light sources to provide an indication of the composition of the traversed 
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human body tissue. 

18. A method in accordance with Claim 17 wherein the indication of the 
composition of the traversed human body tissue takes the form of a photon absorption 
coefficient spectrum for the traversed human body tissue. 

19. A method in accordance with Claim 17 wherein the light photons from the 
different light sources are applied in a sequential manner. 

20. A method in accordance with Claim 17 wherein the human body tissue is 
female human breast tissue, the delivery and collection points are located on the 
surface of a female human breast and the resulting indication of tissue composition can 
be used to detect an abnormal tissue condition. 



31 



WO 03/012391 



PCT/US02/23628 



1/9 




FIG. 1 



WO 03/012391 



PCT/US02/23628 



2/9 



FDPM 



26 



jli^ (X)=AX 



■28 



FIG. 2 



-24 



SS 



ol-RQC) 



S-25 



\y~ 27 



■29 



X10~ 



i 



Q3 




750 800 850 

Wavelength(nm) 
FIG. 3 



900 950 



WO 03/012391 



PCT/US02/23628 



3/9 




U > i i i i i i 

650 700 750 800 850 900 950 1000 

Wavelength(nm) 
FIG. 4 




OA 1 1 . 1 , 1 

650 700 750 800 850 900 950 1000 

Wavelength(nm) 
FIG. 5 




650 700 750 800 850 900 950 1000 



Wavelength(nm) 
FIG. 6 



WO 03/012391 



PCT/US02/23628 



4/9 



1.15 i 



1.1- 




0 75 ~ 'I.,, 

650 700 750 800 850 900 950 



Wavelength(nm) 
FIG. 7 



WO 03/012391 



PCT/US02/23628 



5/9 






WO 03/012391 



PCT/US02/23628 





Wavelength (nm) 
FIG. 12 



WO 03/012391 



PCT/US02/23628 



7/9 




650 700 750 800 850 900 950 1000 

Wavelength (nm) 
FIG. 13 



I 



03 



0.018 -f 
0.076 
0.014 - 
0.012 - 
0.01 - 
0.008 
0.006 



0.004 



650 700 750 800 850 900 950 1000 

Wavelength (nm) 
FIG. 14 



WO 03/012391 



PCT/US02/23628 



8/9 




700 800 " 900 1000 

X(nm) 
FIG. 16 



WO 03/012391 



PCT/US02/23628 



9/9 




800 900 

Wavelength (nm) 
FIG. 17 



1000 



X10~ 



I 



CO 




700 



800 900 

Wavelength (nm) 
FIG. 18 



1000 



